Surface Superconductivity of Dirty Two-Band Superconductors: Applications to MgB 2 . 
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The minimal magnetic field H C 2 destroying superconductivity in the bulk of a superconductor is 
smaller than the magnetic field H c3 needed to destroy surface superconductivity if the surface of 
the superconductor coincides with one of the crystallographic planes and is parallel to the external 
magnetic field. While for a dirty single-band superconductor the ratio of H c3 to H C 2 is a universal 
temperature- independent constant 1.6946, for dirty two-band superconductors this is not the case. I 
show that in the latter case the interaction of the two bands leads to a novel scenario with the ratio 
Hc3/H C 2 varying with temperature and taking values larger and smaller than 1.6946. The results are 
applied to MgB 2 and are in agreement with recent experiments [A. Rydh et al., cond- mat/0307445]. 

PACS numbers: 74.20.De, 74.25. Op, 74.81.-g 
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Introduction. It is well-established that strong mag- 
netic field destroys superconductivity. If an external field 
H applied to a type-II superconductor exceeds the sec- 
ond critical field H C 2, the bulk order parameter in the 
superconductor vanishes. However, even for H > H c % 
superconductivity might still exist in a thin layer close 
to the surface if H is smaller than the third critical field 
H c3 > H c2 [l]. In this paper I investigate the onset of 
superconductivity via surface nucleation for the field H 
slightly below the threshold H c ^. 

In their pioneering work[l] Saint- James and de Gennes 
have shown that if the external magnetic field is ap- 
plied parallel to the surface of an isotropic single-band 
superconductor [2] with a temperature close to the tran- 
sition temperature T c , the ratio UcilUci takes the uni- 
versal value r\ = 1.6946 independently of the supercon- 
ducting material. It turns out[l] that for H slightly 
below H C 3 the superconducting order parameter exists 
within the distance C(T) (the coherence length of the 
superconductor) from the surface. For distances exceed- 
ing £(T) the order parameter approaches zero rapidly. 
The dependence of the ratio H C i/H C 2 on the material 
properties [3, 4, 5], sample geometry and topology [6, 7, 8], 
and temperature [9, 10, 11] has become a subject of in- 
tensive investigations. 

A novel window for investigating surface supercon- 
ductivity was opened after the discovery of the two- 
band superconductor MgB 2 [12]. Not only has it a rel- 
atively high (ss 40 K) T c but also there exist two dif- 
ferent superconducting gaps. As the consequence of this 
fact, various properties of MgB 2 are quite different from 
those of single-band superconductors. For example, the 
anisotropy 7 (T) = H^{T)/H^{T) (here, H { c f(T) 

(c) 

and H c2 (T) stand for the second critical fields in the 
ab and c-directions respectively; note that the crystal of 
MgB 2 is uniaxial) of the second critical field H C 2 exhibits 
strong dependence on temperature, see e.g. Ref. [13]. For 
single-band superconductors this ratio is constant. 

Another puzzle is that j(T) varies widely in different 
experiments[14]. This can be attributed to the existence 




FIG. 1: Hcz/Hc2 for the surface of the MgB 2 crystal co- 
inciding with the a&-plane and the external magnetic field 
lying in the afe-plane for different ratios of the diffusivi- 
ties D2,c/D\, c — 100,300,600 in the two bands (dash-dot, 
solid, and dashed respectively) as a function of T/T c . Inset: 
H c3 /H c2 close to T/T c = 1. For T = T c , H c3 /H c2 = 1.6946. 



of surface superconductivity which might affect the ob- 
servable values of H c2 and, hence, the anisotropy. Conse- 
quently, the determination of the third critical field H c ^ 
is a very important problem. In a recent experiment [15] 
it has been shown that H C ^/H C 2 for MgB 2 might be re- 
duced. 

In the present paper I investigate the ratio HcijHci 
for a dirty MgB 2 crystal. The existence of two different 
gaps manifests itself through the remarkable dependence 
of H c s/H C 2 on temperature. This is in sharp contrast 
with the case of a dirty single-gap superconductor where 
H c3 /H c2 = i] in the whole temperature range. For a 
magnetic field lying in the a6-plane of the MgB 2 crystal 
I find that if one starts decreasing temperature, H C 3/H C 2 
first exhibits a maximum at T « 0.99T C and then a min- 
imum at T ps 0.9T C . As temperature decreases further, 
H c z/H C 2 increases and tends to a value slightly below rj, 
see Fig. 1. Naively, one could try to use the Ginzburg- 
Landau theory (GLT) in order to find H C ^/H C 2. However, 
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as it will be explained below, this would lead to the ratio 
H c s/H C 2 — rj in the whole temperature range, i.e. one 
needs a more rigorous approach in order to explain the 
deviation of H c3 /H c2 from rj. 

General formalism. An appropriate tool to investigate 
magnetic properties of dirty superconductors is the Us- 
adel equations[16]. For two-band superconductors they 
have been derived by Koshelev and Golubov[17] and by 
Gurevich[18]. Since I investigate the onset of supercon- 
ductivity near iJ c2 or H c3 , it is possible to write the Us- 
adel equations in the linearized form 

A a = 2irT (2) 

Here, uj = 2i:T{n + 1/2), n = 0, 1, . . . and ojd are the 
Matsubara and cutoff phonon frequencies. D a j is the 
diffusion coefficient of the band a = 1, 2 along the direc- 
tion j = x,y,z. The indices 1 and 2 correspond to the 



a- and 7r-bands respectively. A is the vector potential. 
A a and f a are the superconducting gap and anomalous 
green function for the band a. The matrix A represents 
the strength of the coupling parameters and has the val- 
ues An « 0.81, A 22 = 0.285, X 12 ~ 0.119, and A 2 i « 0.09, 
see Rcf. [19]. In the present paper I will concentrate on 
two geometries: i) the magnetic field is parallel to the 
c-axis and the surface of the crystal; ii) the surface of 
the superconductor coincides with the a6-plane and the 
field H lies in the a6-plane. As I will show, in case i 
H c z/H C 2 = v at any temperature. In case ii, H c3 /H C 2 is 
shown in Fig. 1. I assume that H is parallel to the 2-axis. 



Choosing the gauge as A y = Hx I look for the so- 
lution of the form f a = f a (iv,x) exp (ik y y + ik z z) and 
A Q = A a (x) exp [ik v y + ik z z). In general, Eqs. (1) and 
(2) define a sequence of solutions corresponding to dif- 
ferent eigenvalues H = H c2 or H = H c3 . One should 
look for the maximal possible values of H C 2 or H c3 . This 
corresponds to the case k z = 0. Substituting the Ansatz 
for f a and A a into (1) and (2) I obtain the system of 
equations 



A" 



Ai(a;) 
A 2 (x) 



I 2ttT £ 



U1>0 



uj+H 1 (x ) 



2nT E 



Ai(x) 
A 2 (o;) 



(3) 



^0 u+H 2 (x a ) J 



where 

H M = ~^g-2+^[-^) (X-X ) , (4) 

with D a ^ x — D a , z — D a a for case i and D a x = D a a 
and D a z = D a c for case ii (D a a = D a j, ^ D a c 
are the diffusion coefficients along the crystallographic 
axes), xq = fc y $o/27riJ is the parameter characteriz- 
ing how far away the superconducting nucleus is situated 
from the surface. Note that xo is the same for the both 
bands. The operator (4) can be rewritten in the form 

a {x'o), with 

d 2 2 
ha{x' a ,o) — ~ q x , 2 + ( X 'a ~ x 'a,o) ! (5) 

where I have made the variable substitution 

x = (3 a x' a and x = f3 a x' a0 , with j3 a = 

( J D a ,,/^ Q ,,) 1/4 ($o/2^) 1/2 . 

The system (3) should be solved with the bound- 
ary conditions (BC) dA a /dx\ x =o = 0, and A a (x — > 



+00) — > 0, a = 1,2 valid for geometries i and ii, see 
above. For H c3 the BC are well-established for dirty 
superconductors [5]. For H C 2 the application of these BC 
gives the same result as the BC requiring the appearance 
of a superconducting nucleus in the bulk. The procedure 
of finding H c3 /H C 2 is following: first, set xo = in (3) 
and find the maximal possible field H for which the so- 
lution satisfying the BC exists. This gives H c2 . Next, 
for xo =/= find the maximal field H = H(xo) for which 
the solution of (3) exists. Then, H c3 = max Xo {H(xo)}- 
I would like to mention that there are complementary 
approaches for calculating iJ c2 based on macroscopic 
theory [20, 21] and GLT[22, 23]. 



Here, it is instructive to study briefly the case of a 
single-gap superconductor. This corresponds to Ai 2 = 
A 2 i = 0. H C 2 and H c3 are then determined by those for 
band 1 (as An > A 22 ). The solution for A\(x) is propor- 
tional to the ground state wave function of the operator 
Hi(xq) and A 2 {x) = 0. Substituting this Ansatz into (3) 
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I obtain the transcendental equation of the form 



1-AnE 



M) uj + (ttH/Qq) yj Di, x Di, z € Q (x' 1i0 ) 



0, (6) 



with eo(x' a ) the lowest eigenvalue of the operator 
h a (x' a ). The field H C 2 can be found as the solution of 
the above equation for x' l = 0; note that eo(0) = 1. As- 
sume, a certain value of the magnetic field H c2 is found; 
let us change the parameter x\ . This leads to the de- 
crease of the eigenvalue cq(x' 1 )[1]- In order to satisfy 
Eq. (6) one has to increase the field H; that is why 
H c3 > H C 2- The minimal £0(^10) can be realized for 
x' 10 = 0.7618[1] and is equal to 0.5901 [1]. This means 
that H c3 /H c2 = 1/0.5901 = r/ for any temperature T. 
Remarkably, it is not necessary to solve Eq. (6) in order 
to find the ratio H c3 /H C 2 1 although the determination of 
H c2 or H c3 alone would require the complete analysis. 

Case i. In this case, the ratio H c3 / H C 2 = 77 at all tem- 
peratures. This is a consequence of the fact that the op- 
erators Hi(xo) and H 2 (xo) have identical eigenfunctions 
(since Di, x /Di, z = D2, x /D 2 , z = !)• The functions Ai(x) 
and A2(x) are proportional to the ground state wave 
function of the operator Hi(xo) (or #2(^0))- The equa- 
tion determining the critical fields H C 2 and H c3 has the 
form F(TrD ha He (x' lfi )/<S> , ■KD 2 , a He {x 2fi )/<£> ) = 0, 
with F(yi, 7/2) a certain function of two arguments. Note 
that in the present case (3\ = ($2 and, consequently, 
x[ = x' 2 and eo(x' 1 0) = e o( x 2 o)- The maximal value 
of H c3 can be realized for x[ = 0.7618 and is equal to 
i]H C 2 for all temperatures. 

Case ii. If the magnetic field lies in the a6-plane, 
the operators H\(xq) and H 2 (xo) have different eigen- 
functions. This leads to a complicated transcendental 
equation depending on all eigenvalues of the operators 
-Hi(xo) and H 2 (x ) and not only on the ground state 
ones. Eqs. (3) can be solved via expanding functions 
Ai(x) and A2(x) over the eigenfunctions of the opera- 
tors Hi(xq) and ^(xo)- I have truncated the basis of 
the operators H\{xq) and i?2(^o) to subspaces consisting 
of 70 eigenfunctions and solved the system (3) numeri- 
cally. In Fig. 1 I show the results of the numerics. Here, I 
take D lia /D liC = 40.0, L> 2 ,a/-t>2, c = 0.665. These ratios 
can be obtained using the results for the average veloc- 
ity on the MgB 2 Fermi surfaces and assuming isotropic 
scattering, see Refs. [24, 25]. The ratio D2. c /Di :C takes 
three values: 100, 300, and 600. This choice is moti- 
vated by the facts that the ratio £ ) 2,c/£'i,c ~ 100 can 
be obtained assuming that the scattering rate of elec- 
trons is the same in both bands. On the other hand 
R = D 2 , c /Di,c = 600 gives a better fit with experiments 
on the anisotropy measurements [24]. 

The results are shown in Fig. 1. For temperatures 
T < 0.6T C I have found that the ratio H c3 /H c2 is nearly 
constant and has a value slightly below 77. This can 
be explained by the fact that at low temperatures the 



fields H C 2 and H c3 are determined mostly by the a- 
band whose coherence length is much smaller than that 
of the 7r-band. At low T the magnetic field H c2 de- 
pends on the ground state eigenvalues of the operators 
H\{xq) and H2(xq) and the contribution of excited states 

(~ - \ 1/4 
is negligible[17, 18]. The ratios [Di yX /Di yZ ) and 

/- \ 1/4 

D2,x/D2 t z) determining the length xo are equal to 
w 2.5 and 0.9 respectively. This means that one can 
maximize the field H c3 by choosing x' 10 = 0.7618. The 
length x' 2 then is large and the ground-state eigenvalue 

of the operator ^(a^.o) ^ s cl° se to 1. The ratio H c3 /H c2 
then can be calculated as follows: take the zero tem- 
perature expression for iJ C 2[18, 24] and make there a 
substitution Dij — > D\^jr\. At low temperatures [18, 

24], H c2 = $oT c cxp(.g/2)/2 7 (C 1:Q C liC i?2,a^2, c ) 1/4 , 
with g = (Xl/w 2 + ln 2 k/4 + 2A_ In k/w) 1/2 - X /w, 

K = D 2,a-D2,c/-Dl,o-Dl,ci A_ = An — A22, Aq = 

1 /2 

(X 2 _ + 4Ai 2 A 2 i) , w = A11A22 - A12A21, and In 7 = 

0. 5772. This procedure yields H c3 /H c2 = 1.688,1.690, 
and 1.691 for D 2 , c /D hc = 100, 300, and 600 respectively. 
The values obtained in the numerics are slightly larger 
(but still smaller than 77) due to a small contribution to 
the adjustment of xo from the 7r-band. 

If one increases temperature, the ratio H c3 /H C 2 de- 
creases and exhibits a minimum at T ~ 0.9T C . Then, 
the value of H c3 / H c2 goes up and takes a maximum at 
T ~ 0.99T C . At T = T c , H c3 /H c2 = r). The nontriv- 
ial behavior of the ratio H c3 /H c2 in MgB 2 is due to the 
changing relative importance of the 7r-band. While at 
low T it is unimportant, at high T it gives a comparable 
with the tr-band contribution to the fields H C 2 and H c3 . 

Let us analyze the situation close to T c in more detail. 
In particular, let us explain why at T = T c , H c3 /H C 2 = r\. 
For T c — T <C T c the field H c3 is small and so are 
the eigenvalues of the operators Hi(x ) and H2(xo), 

1. e. one can use the expansion J2Z>o + H a (x )j w 

E">o l/w - E^>o H a (x )/u; 2 + ... and Eqs. (3) can be 
rewritten in the form 



W 



A 1 (x) 
A 2 (x) 



Ri(x ) 






#2(2:0) 



Ai(!B) 

A 2 (x) 



(7) 



with W 



A" 



- ln(27^ D /7rT) 1 2 and R a (x ) = 
2nT^ l/(oj + H a (xo))-2nT^>o^- T c is dc- 
termined by det W = 0. Solving the system (7) for Ai(x) 
I obtain 

(w n R 2 {xo) + W 22 Ri{xo) - R2(*o)Ri(xoj) Ai(ar) = 0. 

(8) 

Eq. (8) allows to determine the fields H c2 and H c3 in a 
regular way. To lowest order, one can neglect the term 
R 2 (x )Ri (x ). 
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The equation (WuR 2 (x ) + W 2 2Ri(x f) Ai(x) = 

has the ground-state solution of the same form as Eq. (3) 
for a single-gap superconductor (the case Ai 2 = A21 = 0) 
with D 1>x -> D x = WpDi tX + Wnb 2 , x and D 1>x -» 
Dz = W 22 Di z + WuD 2z and the problem of finding 
H C 3 becomes equivalent to the original one considered 
by Saint- James and de Gennes[l]. Consequently, to low- 
est order in T c the ratio H c ^/H c2 has the same value 77 
as in the case of a single-gap superconductor. The ap- 
proximation described above is equivalent to the GLT. 
Consequently, the GLT is unable to explain deviations of 
H c3 /H c2 from r\. 

The maximum of H c3 /H c2 takes place very close 
to T c and is at the boundary of the accuracy of the 
present numerical calculations. Hence, an analytical ap- 
proach would be useful. Temperature corrections to 
H c z/H c2 can be found by expanding operators i?i(xo) 
and R 2 {xq) to second order in (T c — T). One can de- 
compose the operator in the left-hand side of (8) as a 
sum L\(H, xq) +L 2 (H, xq), with Li(H, xq) ctT c — T and 
L 2 (H,x a ) oc (T c — T) 2 . Let \cj) ) be the solution of the 
equation Li(H, Xo)\(t>o) = and H = H^(x a ) the crit- 
ical field to this order. The correction to the eigenvalue 
can be found perturbatively and are determined by the 
implicit relation 

(<h\Li{H,x )\<h) + (<j> \L 2 (H^,x )\(Po) = 0. (9) 



To this order x = 0.7618(Dx/D z ) 1/4 (l/2TrH^) 
Temperature corrections to H c ^ due to change in xq are 
proportional to (T c — T) 3 and can be neglected. Since 
|</>o), Li(H,xq), and L 2 (H,xq) are known, one can find 
H C 3/H c2 analytically. Straightforward but quite cum- 
bersome calculations [26] show that near T c , H C 3/H c2 « 
77 + b(T c - T)/T c , with b ~ 1 for D 2>c /D ltC in the range 
from 100 to 600, in accordance with the numerics, see 
inset in Fig. 1. 

Experiment. Recent experiments [15] show that that 
the ratio H c ^/H c2 is reduced in the case ii. The values 
H c z/H c2 « 1.5 in the temperature range 20-30 K have 
been reported. For the case i H c ^,/H c2 w 1.7[14, 15]. 
The present theory gives that H c ^/H c2 — r\ for case i 
and Hcs/H c2 < r\ for case ii, in agreement with [14, 15]. 
Theoretical calculations show[24] that the anisotropy of 
MgB 2 is distributed over a wider temperature range in 
an experiment that theory suggests. A somewhat simi- 
lar situation takes place in the present work for the ratio 
H c z/H c2 . There are two main sources of deviation be- 
tween theory and experiment. First, surface quality[27] 
might affect H c3 /H c2 . Second, MgB 2 is situated some- 
where at the boundary of the applicability of the weak- 
coupling BCS-theory. It would be very interesting to 
repeat the calculation done in the present paper starting 
from the Eliashberg equations. 

The method described above can be generalized for 



,1/2 



an arbitrary direction of crystallographic axes with re- 
spect to the surface of a superconductor. For strongly 
anisotropic superconductors, surface superconductivity 
might disappear if the surface does not coincide with 
crystallographic planes [8]. Estimates [26] show that 
MgB 2 is sufficiently anisotropic in order to observe this 
kind of effects. The detailed analysis of the onset of sur- 
face superconductivity in this case is challenging for both 
theorists and experimentalists. 

In conclusion, I have presented the calculation of the 
ratio H c3 /H c2 for the two-band superconductor MgB 2 
in the dirty limit. Remarkably, in contrast to the 
case of a single-gap superconductor, the above ratio is 
temperature-dependent. The Ginzburg-Landau theory 
is unable to explain deviations of H c3 /H c2 from r\. 

The present work is supported by the Packard Foun- 
dation. D. A. G. thanks M. Angst and A. E. Koshelev 
for helpful discussions. 
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